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Abstract 

In the first part of the talk, I review the applications of loop equations to the 
matrix models and to 2-dimensional quantum gravity which is defined as their con- 
tinuum limit. The results concerning multi-loop correlators for low genera and the 
Virasoro invariance are discussed. 

The second part is devoted to the Kontsevich matrix model which is equivalent 
to 2-dimensional topological gravity. I review the Schwinger-Dyson equations for 
the Kontsevich model as well as their explicit solution in genus zero. The relation 
between the Kontsevich model and the continuum limit of the hermitean one-matrix 
model is discussed. 
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1 Introduction 



The relevance of matrix models to the problem of genus expansion of Feynman graphs 
goes back to the original work by 't Hooft ||Hoo74 1 . An explicit solution for the simplest 
case of the hermitean one-matrix model had been first obtained by Brezin, Itzykson, Parisi 
and Zuber [BIPZ78] in genus zero and then extended in [ Bes79| , IZ80 , BIPZ80 to next 
few genera. 

The modern interest in matrix models is associated with the context of statistical theo- 
ries on random lattices and discretized random surfaces [|Kai85 , |Dav85| , |ADF85| , [KKM85 



as well as with the conformal field theory approach to 2D quantum gravity [po!871 , [KPZ88 



Dav88|, |DK89|| . A connection between continuum limits of the matrix model and minimal 



conformal models had been conjectured by Kazakov ||Kaz89|1 on the basis of genus-zero 
results. 



The whole genus expansion of 2D quantum gravity has been constructed in [ BK90 
DS90| , |GM90b|| taking the 'double scaling limit' of the (hermitean) one-matrix model. 
Moreover, the specific heat turns out to obey a (non-perturbative) equation of the 
Korteweg-de Vries type so that a relation between the continuum limit of the matrix 
models and integrable theories emerges ||GM90a| , BDSS90| . 

While these results were obtained using orthogonal polynomial technique, one more 
method — that of loop (or Schwinger-Dyson) equations — is custom in studies of ma- 
trix models. Loop equations had been proposed originally for Yang-Mills theory both 
on a lattice |[Foe79| , pgu79| , [Wei79|| and in the continuum ||MM79| , |Pol80|| (for a review, 
see |[Mig83|| ) and then were applied [pR80| , [Fri81| , [Wad8 1| to matrix models. A modern 
approach to loop equation which is based on its interpretation as a Laplace equation on 
the loop space can be found in ||Mak88| , |Mak89| , [HM89|1 . The recent applications of loop 
equations to 2D quantum gravity have been initiated by Kazakov [ |Kaz89|| . The role of 
loops in 2D quantum gravity is played by boundaries of a 2-dimensional surface. 

Ref . |[Kaz89| deals with genus zero. The whole set of loop equations for 2D quan- 
tum gravity was first obtained by David 



Dav90 



taking the 'double scaling limit' 
of the corresponding equations for the hermitean matrix model. As was shown in 
Mig83| , Pav90| , |AM90|| , these equations can be unambiguously solved order by order 
of genus expansion. However, this solution is non-perturbatively unstable | pav9(J| ] as it 
should be for 2D euclidean quantum gravity. 

One of the most interesting results which are obtained with the aid of loop equations 
is the fact that the partition function of 2D quantum gravity in an external background 
is the r-function of KdV hierarchy which is subject to additional Virasoro constraints 
J91| , DV V91 1 . This proves a conjecture of Douglas |pou90 |. The existence of Virasoro 



algebra was extended to the case of the matrix model at finite N in ||AJM90|, |GM1VP9T 



IM91| , |MM90|| while the relation to the continuum Virasoro algebra of |[FKN91| , PVV91 
had been studied by M 4 |MMMM91|| . 

The second application of loop equations concerns the relation between 2D quantum 
and topological [|LPW88| , |MS89|1 gravities. As Witten [|Wit90|l had conjectured, these two 
gravities are equivalent. This conjecture has been verified in genus zero and genus one 



for 2D quantum gravity coincide with the recursion relations between correlators in 2D 



topological gravity which were obtained by Verlindes |[VV91 



From the mathematical point of view, a solution of 2D topological gravity is equivalent 
Wit90|| to calculations of intersection indices on the moduli space, M g>s , of curves of genus 



g with s punctures. Interesting results for this problem have been obtained recently by 
Kontsevich |[Kon91|| who has represented the partition function of 2D topological gravity 
as that of a (hermitean) matrix model in an external field. It is worth noting that 
the Kontsevich matrix model is associated with the continuum theory. Therefore, it 
should be directly related to the 'double scaling limit' of the standard one-matrix model 
PR9q , PB9q , rGM90b|| . 

The Kontsevich model can be studied by the method of loop equations. As has been 
shown recently by Semenoff and the author [ |MS91a|| , the Schwinger-Dyson equations for 
the hermitean one-matrix model in an external field, which is equivalent to the Kontsevich 
model, can be represented as a set of Virasoro constraints imposed on the partition func- 
tion. The large- N solution of these equations, which is known from the work of Kazakov 
and Kostov [ |KK89| | , solves the Kontsevich model in genus zero | |MS91a| | showing explicitly 
the equivalence of 2D topological and quantum gravities to this order. 

The fact that the partition function of the Kontsevich model obeys the same set 
of Virasoro constraints |[FKN91| , |DVV91|1 as the continuum limit of the hermitean one- 
matrix model has been proven recently by Witten [ Wit9l 1 using diagrammatic expansion 
and by (A.M.) 3 [MMM91|1 using the Schwinger-Dyson equations. This demonstrates an 
equivalence of 2D topological and quantum gravities to any order of genus expansion. 

In the first part of the talk, I review some works ||AM9Q , |AJM90| , |Mak90| , |MMMM91|j on 
applications of loop equations both to Ax N matrix models at finite N and to 2D quantum 
gravity which is defined as their continuum limit. The results for multi-loop correlators 
in low genera and the Virasoro invariance both at finite N and in the continuum are 
discussed. The second part is devoted to the Kontsevich matrix model which is equivalent 
to 2-dimensional topological gravity. The Schwinger-Dyson equations for the Kontsevich 
model as well as their explicit solution in genus zero [ |i\lS91a|| is reviewed. Some original 
results concerning the Kontsevich matrix model are reported. 



2 Matrix Models and 2D Quantum Gravity 



2.1 Loop equation for hermitean matrix model 

The hermitean matrix model is defined by the partition function 

= J VM exp-trV(M) (2.1) 
where M is the N x N hermitean matrix. V stands for a generic potential 

oo 

V{p) = Y. t *P k - ( 2 - 2 ) 

k=0 

The coupling plays here the role of a source for the operator trM fc while V(p) is a 
source for the Laplace image of the Wilson loop tr [l/(p — M)] : 

/+2OO+0 rjrn 1 
— V(p) tr -. (2.3) 
-ioo+o 2ni KFJ p-M K J 

The correlators (trM fcl ... tiM km ) c , where the average is defined with the same 
measure as in Q2.1|), can be obtained differentiating \ogZ^ w.r.t. t kl , . . . , t fcm while loop 
correlators can be obtained by applying^ 

r) 00 8 



sv(p) t f„ at k 

so that the m-loop correlator reads 

1 1 6 6 

W u ( Pl , ...,p m ) = (ti . . . tr ) c = —— . . . —— log Z». (2.5) 

Pi - M p m - M oV^px) 5V{p m ) 

To calculate the actual values for the given model, say for the matrix model with cubic 
potential, one should, after differentiations, put t^s equal their actual values, say tj. = 
for k > 3 in the case of cubic potential. 

The loop equation can be derived using the invariance of the integral under an (in- 
finitesimal) shift of M and reads 

/ ^P^W^u) = (W^p)f + -^W^p). (2.6) 
JCi 2m (p — to) oV [p) 

The contour C\ encircles singularities of W u (uj) so that the integration is a projector 
picking up negative powers of p. Eq.(2.6) is supplemented with the asymptotic condition 

pW n {p) — > N as p -> oo (2.7) 

which is a consequence of the definition ( |2.5|) . 

Notice that one obtains the single (functional) equation for W u (p) . This is due to the 
fact that tr V(M) contains a complete set of operators. Such an approach is advocated in 
DVV9 1 1 [Mak90|l . The set of equations for multi-loop correlators (|2.5| ) , which is considered 



in 



Dav90| , FKN91 |, can be obtained from Eq.( |2.6| ) by m — 1 -fold application of 5/5V(pi 



The system of the standard Schwinger- Dyson equations for the connected correlators 
(trM fcl . . . trM km ) c can be then obtained by expanding in powers of p 1 , . . . ,p„^ . 



2.2 Solution in 1/N 



Eq.( |2.6|) can be solved order by order of the expansion in 1/N 2 (the genus expansion). 
The second term on the r.h.s. represents the connected correlator of two Wilson loops 
and is, in our normalization, of order 1 as N — > oo while two other terms are of order 
N 2 since W n (p) and V(p) are of order N . Therefore one can omit it as N —>■ oo (which 
corresponds to genus zero = the planar limit). 

The simplest (one-cut) solution of Eq.( [2.6|) as N — » oo reads [ Mig83| 



where 



2W H{0) (p) = V'(p) - M(p)yf(uj-x)(u-y) 
du V'{p) - V'(u) 



M(p) 



Ci 2m (j) _ oj)JJoj — x){uj — y) 



(2.1 



(2.9) 



is a polynomial of degree K — 2 if V{p) is that of degree K. The ends of the cut, x and y, 
are determined from the asymptotics ( |2.7p : 



duo 



& 2m J(oj-x)(u-y) 



2N 



du 



uV'(u) 



Ci 2m /( LJ _ x )( LJ _ y ) 



W(x,y). (2.10) 



For the even potential (V(— p) = V(p)), the first of these equations yields y = —x = y/z 
which simplifies formulas. This case is called the reduced hermitean matrix model. 

The multi-loop correlators in the planar (genus zero) limit can be obtained by varying 
according to the r.h.s. of Eq. ( |2.5| ). The 2-loop correlator reads | |AJM9U[ 



W H{0 \p,q) 



1 



2pq — (p + q)(x + y) + 2xy 



A (P Q) 2 [yj(p-x)(p- y)yj (q -x)(q- y) 



(2.11) 



while an expression for the 3-loop correlator is given in [|AJM90|| . Note that the 2-loop 
correlator (|2.11[) depends on the potential, V, only via x and y but not explicitly. This is 
not the case for all other multi-loop correlators. 

To calculate 1/N 2 correction to ( |2.8|) one needs W n ^°\p, p) which enters the r.h.s. of 
Eq.(U). Eq.(plD yields 



W H(0 \p,p) 



(x - yf 



16(p — x) 2 (p — y) 2 
and one can now obtain W/ 11 ^) (p) by an iteration of Eq. ( |2.6| ) 



(2.12) 



du 



The result 

(x - y) 2 



yj(p-x)(p-y) 2ni ( u - P) M (u) 16(o; - x) 2 (uj - yf 



(2.13) 



is unambiguous ||Mig83| , |Dav90|| provided that one requires analyticity of W m ^ l \p) at 
zeros of M(p). This procedure of iterative solution can be pursued order by order of 



2.3 Continuum loop equation 

The continuum limit of the reduced hermitean matrix model is reached as N — > oo while 
K-l conditions W {n) (z c ) = (W(z) = W(—s/z, xfz)) with n = 1, . . . , K-l are imposed 
on the couplings, tk, in addition to ( |2.10| ) at K th multi-critical point. 2D quantum gravity 
corresponds to K = 2. The 'double scaling limit' can be obtained if one expands around 
the critical point: 

, an , ay/X 

V -> ^7= I * -> V^c- V7 =, (2.14) 

so that 7r and A play the role of continuum momentum and cosmological constant, re- 
spectively. The dimensionful cutoff a should depend on N such that the string coupling 
constant G = jV^a" 2 *- 1 would remain finite as N -> oo lEKSn, IDSUUl CMOT. 



To obtain the continuum limit of loop correlators fl2.5|) , it is convenient to introduce 
the even parts 

■-«") s W^Hn) ' ' ' ' 0g Z """ (2 ' 15) 

where Z r ^ and V even (p) means, respectively, (|2.1| ) and (|2.2|) with t2fc+i = : 

W even (pi, . . . ,p m ) differs from W (j>x, ■ ■ ■ ,p m ) by correlators of products of traces of odd 
powers of M. Near the critical point, one gets 

W H (pi, ...,p m )^ 2 m - 1 W even (p 1 , . . . ,p m ). (2.17) 

This formula can be proven analyzing loop equations or by a direct inspection of multi- 
loop correlators ||AJ M9U| ] . 

The continuum loop correlators can be obtained by the multiplicative renormalization 

Pa^gq , imgq , groggg 

W*(pi, . . . - 2 m a- m G^ m ~ 1 W cont (7i 1 , ...,n m ) for m > 3 (2.18) 
while additional subtractions of genus zero terms are needed for m = 1 and m = 2 : 

Wgr(p) - \v\p) -> -±=(2W cont (n) - J'(tt)) , (2.19) 
2 avG 



W^v(pi,P2) -> 4a- 2 iy co „i(vr 1 , vr 2 ) + / . (2.20) 

^ (V^l + V 7r 2) V 71 "!^ 

For latter convenience, W^(pi, . . . ,p m ) on the l.h.s.'s of these formulas is the multi-loop 
correlator for the 2N x 2N reduced hermitean matrix model. 

J {it) on the r.h.s. of Eq.( |2.19) ) plays the role of a source for the continuum Wilson 
loop: 



oo r oo q 

J(tt) = V T n7 r n +5 , = -Y vr-"- 3 / 2 -^- (2.22) 

with Tfc being sources for operators with definite scale dimension. Therefore, Eqs. fl2.18p , 
( |2.20| ) can be derived from Eq. fl2.19|) by varying w.r.t. J(n). 

The continuum loop equation can be obtained from fl2.6|) substituting fl2.19|) , fl2.20|) : 

r dQ J' (to) 2 8W cont (*) G T 2 

/ T~7 fj7 Ww(«J = (^Ww) + G + — -J + — . (2.23) 

JCi 27U (7T — \l) oJyTT) lO/T^ 16/T 

This equation describes what is called the 'general massive model'. It corresponds to 
arbitrary J(n) and interpolates between different multi-critical points. For K th multi- 
critical point, one puts, after varying w.r.t. J(vr), T n = except for n = and n = K. 



2.4 Genus expansion 

Eq. fl2.23p can be solved order by order in G (genus expansion) analogously to that of 
Sect. 2. 2 . If J(vr) is a polynomial (T n = for n > K), K — 1 lower coefficients of 
the asymptotic expansion of W C(m t{'^) are not fixed while solving in 1/tt and should be 
determined by requiring the one-cut analytic structure in ir. The continuum analog of 
( gg) reads 

2W2.00 = / ^L/'Wf^ (2.24) 
contK ' Jd 2m (tt - fi) v/nr^ 

where u versus {T} is determined from the asymptotic behavior. 

This asymptotic relation can be obtained comparing 1/ir terms in Eq.( [2.23j ). Denoting 
the derivative w.r.t. x = — T /4 by D , it is convenient to represent this relation as 

2x= [ ^J'{to)DW cont {to). (2.25) 
JCi 2m 

For the ansatz ( |2.24j ), one gets 

DwHSti*) = -7=== ~ -7=- ( 2 ' 26 ) 

Eq. fl2.26 ) can be extended to any genera using the representation ||GM90ai |BDSS90|| 

/ 1 \ _1 1 00 
DW cont (n) = 2(x\ U + u(x)--GD 2 ) \x) ■= = 2 £ R n [u] = 2R(n) (2.27) 



4 ; 1 ' n=1 



where the diagonal resolvent of Sturm-Liouville operator is expressed via the Gelfand-Dikii 
differential polynomials ||GD75|| 

,-n-l ^ r>2 „. n n-l„ 



i? n [u] = 2~ ( g"^ — u — D — D~ uDj -1. (2.28) 

Substituting the r.h.s. of Eq. (E£27l) into Eq. (|2~25p , one obtains the string equation 
PM9Q4 PTJ559(| 



The fact that the ansatz ( 2.27| ) does satisfy Eq. fl2.23| ) is shown in 
aid, one applies the operator 



PW9T 



To this 



A = -^D 4 + (u + n)D 2 + l(Du)D, (2.30) 
lb 2 



which annihilates Wcont^) given by Eq. fl2.27|) , to Eq. fl2.23|) . The result vanishes provided 
u satisfies Eq.( [2.29|) and 

-2— 6 —u = D 2 W cont (n) (2.31) 
oJ(n) 

whose expansion in 1/-7T reproduces the KdV hierarchy du/dT n = DR n+1 [u) . 
Comparing (|2.21| ) and ( j2.31| ), one concludes that 

Z cont = exp |~ j* dy{x - y)u{y) + $(T 1; T 2 , . . .)} (2.32) 

where the perturbative solution of Eq.( ^.29j ), which satisfies u(0) = 0, is chosen and 
an integration 'constant' $(Ti,T 2 , ...) depends on all T's except for T . At the given 
multi-critical point, this $ is unessential so that the r.h.s. of Eq.( 2.32f) coincides with 



the continuum partition function which was obtained in pK90| , |DS90| , |GM90b|| using the 
method of orthogonal polynomials. 

The general procedure of solving Eq. ([2.23| ) order by order in G can be now formulated 
as follows. One should first solve Eq. (|2.29| ) to find u versus x and {T n } (this is perturba- 
tively unambiguous). Then Eq. (|2.27| ) determines DW cont {^) while W cont (ii) itself can be 
can be obtained by integrating 

W cont (n) = 2 [ X dxR(n) (2.33) 



where the integration 'constant' can be expressed via $ entering Eq.( p.32| ). This constant 
becomes unessential for K th multi-critical point when T n = except for n = and n = K 
so that u depends only on the cosmological constant A: 



/2 = xA^jKlf 
(2K + 1)(2K)\T K 

2.5 Multi-loop correlators in 2D quantum gravity 

A formula which is similar to Eq. ( |2.33|) exists for the multi-loop correlators: 



;2.34) 



r x 5 5 

Wconti-Kl, ...,7l m ) =2 dx— r • • • — T-R(vT m ), (2.35) 

J oJ{TTi) f)J(7T m _i) 

where the integration 'constant' depends again on T 1; T 2 , .. . . 

Since R(ir) depends on T's only implicitly via u, the following chain rule can be used 
for calculations 



with 8/8u(x) being the standard variational derivative. The expansion of Eq. fl2.36| ) in 1/tt 
reproduces the standard (commuting) KdV flows. An advantage of Eq.( 2.36| ) is that it 
allows to obtain results without solving the string equation fl2.29| ). Therefore, to calculate 
the multi-loop correlator for a given multi-critical point, one can substitutes the solution 
of Eq.( p.29| ) only for this multi-critical point and should not solve it for arbitrary T's. 

An alternative way of calculating correlators in 2D quantum gravity is to take 
the continuum limit ( ^.14j ) of formulas of Sect. 2. 2 with the aid of the renormalization 
( |2.18|) -( [2.20| ). For the case of pure gravity (the K = 2 critical point), the explicit form of 
W con t(n) is known for genus zero [pav90|1 and genus one |AM90fl : 



5 -T 2 
4 2 



(tt - -VX)\/tt + y/X 



G (7T + f y/A) 

90T 2 A(vr + v/A) 5 / 2 



(2.37) 



where [...]_ means subtraction of the term which diverges as tt — > oo as well as that of 
order 0{tt~^). 

The multi-loop correlators are known ||AM90| , |AJM90|| for genus zero: 



7T! + VK- V 7T 2 + VX 



4(tti - tt 2 ) 2 yjm + Va^7t 2 + 4 ( v 7 ^ + V^) 2 V^ 



(2.38) 



d 



,m—3 



1 



' 11 ^ dA m ~ 3 



,n- 

A i=l (7Tj 



%/A) 3 / 2 , 



for m > 3. 



(2.39) 



As is mentioned above, the additional subtraction is needed only for m < 2. 

Analogous formulas can be obtained for higher multi-critical points. For K 
Eq. ( |2.38| ) remain unchanged while the analog of Eq. ( [2.371) reads |[AM90| 



W. 



(o+i) 



cont 



7 -T 3 
4 3 



[n 2 - -nVX + -A)v tt + 
2 8 



4G 



315T 3 A 3 /2(7r + VA) 5 /2" 



;2.40) 



The above expressions for multi-loop correlators agree with those obtained recently 
PE9T1 , |MMS9lj |MSS9l| , |MS9lb|| for the Liouville theory. 



2.6 Complex matrix model 

The complex matrix model is defined by the partition function 

Z% = J DM'DM 1 exp - tr V even {MM^) (2.41) 



where the integral goes over N x N complex matrices and V even is given by (|2.16 ). As is 



seen from this formula, the complex matrix model resembles the reduced hermitean one. 
However, it differs by combinatorics as well as by the fact that averages of traces of odd 
powers of M do not appear. The model Q2.41D has been studied in [|Mor91| , |AMP91|1 using 



The variables t2k in ( 2.41|) play the role of sources for the operators tr (MM^) k while 
yeven - g a source f or the Wilson loop tr [p/(p 2 — MM^)] (comp. ( [OP ) : 



tr ^"(MMt) 



+IOO+0 



dp 

-o 2-R-i 



The analog of (|2.5| ) reads 

^ c ( Pl ,...,p m ) 



y**»(p) tr 



5 



P 



(p 2 - MMt) 



log ^ 



which leads to the following loop equation for the complex matrix model [ |Mak90 | 



-W^(p). 



(2.42) 



(2.43) 



(2.44) 



Id A-xi (p-oj) 5V even (p) 
This equation should be supplemented with the asymptotic condition same as (|2 . 7| ) . 
Comparison of Eq. ([2.44| ) and loop equation for the reduced hermitean model yields 

1, 



for genus zero: 



W#(Pi 



,Pr, 



,Pr, 



(2.45) 



where the correlators W even for the reduced hermitean model are defined is Sect. 2. 3 . 
Notice that the correlator ior N x N complex matrix model enters the l.h.s. while that 
for 2N x 2N reduced hermitean one enters the r.h.s.. This guarantees the asymptotic 
condition ( |2.7|) . The coefficient 1/2 in Eq. (|2.45|) leads to the following relation between 
the partition functions for genus zero: 



Z% oc 



yreduced 
^2N 



;2.46) 



Due to the relation ( 2.45 ), AW c (p) to the leading order in 1/N is given by the r.h.s. 
of Eq. ( |2.8|) with V replaced by V even and y = —x = \fz . The multi-loop correlator to 
leading order in 1/N can be then calculated by varying according to ( |2.43|) . The analog 
of (PTTQ , (|2~lp reads 



W c ^(p,q) 



2p 2 q 2 — zp 2 



zq 



4(p 2 — q 



2\2 



y/p 2 - z^/q 2 - z 



2pq\ , W c(0 \p,p) 



16p 2 (p 2 



(2.47) 

Moreover, an explicit expression for arbitrary multi-loop correlators exists [|AJM90|| for 
the complex matrix model even far from the critical point (comp.( p.39| )): 

m-3 „ 



W c ^( Pl ,...,p m ) 







W(z) dz J 2zW{z) fJi 2{pf - zfl 2 



n 



for m > 3. 



(2.48) 



As is proven in ||Mor91| , [AMP91| , |Mak90| , [AJM90|| , the complex and hermitean matrix 
models belong to the same universality class in the 'double scaling limit'. This implies, 
in particular, that the continuum limit of all multi-loop correlators coincide with those 
of Sects.2.4,2.5 . They do not coincide, generally speaking, for higher genera far from 
the critical points. Eqs.( [2.45]) , ( [2.46 ) remain valid, however, to arbitrary order of genus 
expansion near the critical points. Using ( |2.17|) , one concludes that W$(pi, . . . ,p m ) has 



2.7 Loop equations as Virasoro constraints 



The loop equation (|2.6|) can be represented as a set of Virasoro constraints imposed on 
the partition function. Eq. (|2.6|) can be rewritten, using the definitions ( ^2|) and ( |2.4|) , as 



1 

Z N n 



E 



where the operators^ 



H 



p 



d 



rH 7 H 







dt k 



k=0 ^ L k+n 

satisfy ||AJM9q , |GM1VT91| Virasoro algebra 



d 2 

E — - — 

0<fc<n 9tk0t n -k 



(n — m)L 



H 

n+m' 



Therefore, Eq.( |2.6| ) is represented as the Virasoro constraints 

I™Z™ = for n> -1. 



(2.49) 



(2.50) 



(2.51) 



;2.52) 



These constraints manifest the invariance of the integral on the r.h.s. of fl2.1|) under the 
shift of integration variable 5M = e ■ M n+1 with n > -1 | |AJM9U| , |MM9U|| . 

It is impossible, however, to make in ( |2.50| ), ( |2.52| ) the reduction to even times. For 
even n, this reduction can be done for the first term on the r.h.s. of ( [2.5U| ) but not for 
the second one. Therefore, there exist no Virasoro constraints imposed on z r ^ duced at 
finite N. 

A set of Virasoro operators built up from the even times, t^k, arises for the complex 
matrix model. The loop equation fl2.44p can be represented as Virasoro constraints 



r C yC 



for n > 0: 



A:=0 



2A- 



8 | y, d 2 

dt2(k+n) Q<k<n ^2kdt2( n -k) 



(2.53) 
(2.54) 



The Virasoro invariance is now related ||AJM9q , |MM90|| to the change SM = e(MM*) n M 
with n > 0. 

Analogously, the continuum loop equation (|2.23| ) can be represented as Virasoro con- 
straints which are imposed on Z cont defined by ( [2.21| ). Using ( |2.22| ), one proves that 
Eq. (|2.23|) is equivalent to the continuum Virasoro constraints ||FKN91| , PVV91 



ncont y 



qr* = 'E(k + i/2)T k 



k=0 



for n > — 1; 

d 2 

frr.f 

0<fc<n- 



9Tk+n n<fc<n-i dTkdT n _h_i 



$0,n 



+ 



1QG ' 



(2.55) 
(2.56) 



The relation between the continuum Virasoro constraints ( |2.55|) , ( |2.56|) and those at 
finite N can be studied |1\U\11\11\IT)T| without referring to loop equations. Introducing 



Ga n+l ' 2 kt 2k T{k + 1/2) 



- 4:NVGa5 n n, 



(2.57) 



or, vice versa, 



kt 



2k 



2N5, 



fe,0 



n>k 



tk—n 



J-V 2 T n T(n + 3/2) 
y/G(n-k)\T(k + 1/2) 



and rescaling the partition function 



yC yG 

Zj n — > Z-^ 



.4, 



\n+m 



r(n + 3/2)r(m + 3/2) 
27r(n + m + l)(n + m + 2)n!m! 



1 _— m— n— 1 



one gets from ( |2.52|) , (|2.53|) 



(2.58) 

(2.59) 
(2.60) 



£C^C _ g for n 



16a" 



for n > 0, 



00 <9 <9 2 

OJ-k+n 0<k<n-l O1 kOl n -k-l 



k=0 



^n,0 ^n,-l^( 

16~ 1QG 



The variables {T} are related to {T} by 



T n = T n + a- 



n 



T n _x - WVGa5 nfi 



(2.61) 
(2.62) 

(2.63) 



n + 1/2 

so that the difference disappears as a — > 0. 

Eqs.( |2.57| ), (p3H) are the standard transition ||Kaz89|| to operators with definite scale 
dimensions in the continuum. The rescaling ( [2.59 ) makes finite as a —>■ so that 
— > Z cont . While the operators £^ tend to C c ° nt defined by Q2.56Q as a — > 0, the 
a — > limit is not permutable with differentiating Z^ w.r.t. T n . This is why C^Z^ are 
nonvanishing (even singular for n > 1) as a — > 0. These terms do not appear RMMMM91 
however, when £*? 's act on 



„— ^ ^ A mrl T m T n / y reduced . y 



cont 



(2.64) 



(comp. ( [2.46D , ( p.59p ). Thus, the l.h.s. of Eq.( ]2~l>4| ) defines the proper continuum partition 
function which is annihilated by the continuum Virasoro operators ( 2.56Q . 



3 Kontsevich Model and 2D Topological Gravity 



3.1 2D topological gravity as Kontsevich model 

The starting point in demonstrating an equivalence between 2D topological gravity and 
the Kontsevich model is the Witten's geometric formulation [|Wit90|1 of 2D topologi- 
cal gravity. In this formulation, one calculates the correlation functions of s operators 
a ni (xi), . . . ,a na (x s ) with definite (non-negative integer) scale dimensions rij, living on a 
2-dimensional Riemann surface £ of genus g. Those are expressed [|Wit90|| via the inter- 
section indices 



where C\{Cu\) is the first Chern class of the line bundle (which is the cotangent space 
to a curve at Xj) over the moduli space, Ai g , s , of curves of genus g with s punctures 
and the integral goes over -M 9)S . The normalization factor YliAf(rii), which is related 
to the normalization of the operators a, is to be fixed below. The r.h.s. of Eq.( [3.i|) is 
non- vanishing only if 

X> i = 3^-3 + s, (3.2) 

i 

i.e. the (complex) dimension of M. g . s . Notice the crucial property of correlators in 
topological theories — those depend only on the dimensions n\, . . . ,n s and genus g but 
not on the metric on £ and, therefore, not on positions of the punctures X\, . . . , x s . 

It is convenient to introduce the set of couplings t n which play the role of sources for 
the operators o n . The genus g contribution to the free energy then reads 



F g [t] = ^exp (j2t n J crj}j 



(3.3) 



while the correlator on the l.h.s. of Eq. (|3.1|) can be obtained by differentiating F g [t] w.r.t. 
t ni , . . . , t Hs since the correlators do not depend on xi, . . . , x s . The total free energy can 
be obtained from the genus expansion 

F[t-\] = ^ 29 - 2 F g [t] (3.4) 

9 

with A 2 being the string coupling constant^. Note, that due to the relation Q3.2p, the 
A-dependence of F can be absorbed by the rescaling of t n : 

F[...,t n ,...;X] =F[...,AK"- 1 )t n ,...;l]. (3.5) 
The Kontsevich approach ||Kon91|| to evaluate F[t; A], given by Eq. ( |3.4|) , is based on 



a combinatorial calculation of the intersection indices on -M 9iS . Let us represent Eq.(|3.3j) 



l 9v 

as 



m = Y. E i/vr ' (3-6) 



S>0«l,.",»s S 



3 We use in this Dart of the talk the Witten's normalization IWit9(J of 2D tonoloeical eravitv (Af(n) = 



where 

F (nu..,n s ) = (a ni ( Xl )---a ns (x s )) g t ni ---t ns . (3.7) 

The last quantity can be interpreted as a contribution from a band graph (or a fat graph 
in Penner's terminology [PenSS]) of genus g with s loops and three bands linked at each 



vertex. These graphs were introduced in quantum field theory by 't Hooft ||Hoo74|| . The 



original Riemann surface with s punctures can be obtained from this band graph by 
shrinking the boundaries of bands (forming loops) into the punctures. 

As is proven by Kontsevich [|Kon9 lfl , 



(3-8) 

where Aj are eigenvalues of a iV x JV hermitean matrix A and the sum goes over the 
connected band graphs with #(vert.) vertices. The product goes over the links of the 
graph. Each link carries two indices i, j which are continuous along the loops while each 
of s traces on the l.h.s. corresponds to the summation over the index along one of s loops. 
The combinatorial factor #(aut.) in the denominator is due to a symmetry of the graph. 

Substituting Eq.flSJl into Eqs.(l3?7|) , (^6|) , identifying 

tr (A" 2 "- 1 ) (3.9) 



X(2n - r 



N{n) 

and making use of Eq. ( |3.5|) , one represents the r.h.s. of Eq. Q3.4|) in the form of the 
logarithm of the partition function 

^W[A;a M JVXe -^ ( 3 - 10 ) 

where the integral goes over the hermitean N x N matrix X. The original normalization 
of |[Kon91|| corresponds to A = — 1. 



There is, however, a subtlety in the identification (|3.9| ). The point is that, for an 
N x N matrix A -1 , tr (A~ fc ) are independent only for 1 < k < N while, say tr (A _iV_1 ) 
is reducible. All tr (A~ 2n_1 ) become independent, as it should be for the sources in 2D 
topological gravity, as N —>■ oo. Therefore 

logZ Kmts [A;X) -*F[t;\]. (3.11) 

only as —>■ oo. The equality (|3.11|) is valid in a sense of an asymptotic expansion at 
large A with each term being finite providing A is positively defined. 

Let us explain the Kontsevich results from the viewpoint of the standard analysis 
of the matrix model ( |3.10| ). Z^ onts [A; A] admits the perturbative expansion in A that 
starts from the term 0(A). This term corresponds to the contribution of three puncture 
operators in genus zero ||Wit90|| : 



+3 



where A 2 in the denominator emerges because of Eq.( |3.9| ). The contribution of a generic 
graph with #(vert.) vertices, #(link) links and s loops is proportional to 



A^AP = (XN) S X 2 ^ 2 (3.13) 

in an agreement with Eq. (|3.8|) . 

Notice that while N — > oo, all terms of the perturbative expansion in A contribute 
to F[t; A] in contrast to the standard large- N expansion by 't Hooft [|Hoo74|| when an 
expansion in emerges so that N = oo corresponds to planar graphs only. The 't Hooft 
case can be reproduced if A ~ AT -1 . Then F ~ A^ 2 while 

W Konts [A; i] = l^\ogZ Konts [A; i] - ^F[t; i] (3.14) 

is finite. Therefore, WKonts[A; 0] = F [t] can be obtained in the 't Hooft planar limit. This 
fact has been utilized in |MS91a|1 to solve the Kontsevich model in genus zero. 



3.2 The Schwinger-Dyson equations 



The Kontsevich model can be studied using the custom methods of solving matrix models. 
Since A in Eq.( |3.10| ) is a matrix, the standard orthogonal polynomial technique can not 
be applied. For this reason, the method of Schwinger-Dyson equations has been applied 
||MS9laL 



IMMM91 



to this problem 

To derive the Schwinger-Dyson equations, it is convenient to make a linear shift of the 
integration variable in the numerator on the r.h.s. of Eq.( |3.1U| ). Modulo an unessential 
constant, one gets 



^Konts [A; A] = n 



AiHiAi + Aj) e~ Z 

i>j 



A 2 



(2A)s 



where 



Z[M] = jvX e tr( " 



x_ 

3 



-MX) 



(3.15) 



(3.16) 



and the Gaussian integral in the denominator has been calculated. 



Z[M] which is defined by Eq. ( p.!6|) is the standard partition function of the hermitean 
one-matrix model in an A^ x N matrix external field M. This external field problem, which 
is analogous to the corresponding problem [BNST , PG80| , PRT81|| for the unitary matrix 
model, has been studied recently in [|GN91| , |MS91a|| . While the representation ( 3.10| ) 
is convenient for constructing the perturbation theory expansion, the partition function 
Z\M\ is convenient for deriving Schwinger-Dyson equations. 

The partition function (|3.16|) depends on A^ invariants, m,, — the eigenvalues of M. 
Let us perform the integral over angular variables in the standard way [ 1Z80| , MchSl 
to express Z[M] as the integral over Xi — the eigenvalues of X. Modulo an irrelevant 
multiplicative constant, the result reads 



where A[m] = rii<j( m « — m j) * s the Vandermonde determinant. 

The Schwinger-Dyson equations result from the following change of variables in ( 3.17Q : 



Xi + e n x? +1 for % = 1, . . . , N and n > — 1 in full analogy to the matrix model without 



external field [|AJM90| , |IM91| , |MM90|| . Noticing that in the integrand can be replaced 
by Tj^: when applied to A[m]Z[m], the set of Schwinger-Dyson equations can be written 
in the' form ||MS91a|| 

L n A[m]Z[m] = for n > — 1 



(3.18) 



with 



E 



'_d_ 

drrij 



+ 



d 



n—k 



dm 



dm ^ 



(3.19) 



It is easy to verify by a direct calculation that these operators obey Virasoro algebra 



in — m)L r 



(3.20) 



The Virasoro generators ( |3.19| ) annihilate the totally antisymmetric function 
A[m]Z[m]. One can easily construct the generators C n which annihilate Z[m] itself. 
Let us introduce for this purpose the 'long' derivatives 



d 

Vj = A-^m]— A[m] 
orrij 







+ E 



i 



(3.21) 



dmt r£ mi - mj 

which commute one with each other. The Virasoro constraints ( |3.18 ), ( |3 . 1 9| ) now take 

C n Z[m] = for n > — 1 



the form |MS91a 
and 



1 



(3.22) 
(3.23) 



^n = E -(v,) n+3 + (v,)" +1 m 2 + - x: E(^) fe (v,r- fe 

i [ fc=0 j'^i 

Due to the commutativity of V's, the order in the last term in unessential. As follows 
from the definition ( |3.21| ), the generators ( |3.23| ) obey the Virasoro algebra commutation 
relations, same as (|3.20|) . 

The Virasoro constraints ( |3.22p , ( |3.23| ) turn out to be equivalent to the following 
equation 

(<9*) 2 + ^ z^-zr idi - dj) - ^ j Z[m] = (3.24) 



L . m; — mj 



which is called in [ |MS91a|| the 'master equation'. As is shown in ||GN91] , [MMM91 
Eq.( 3.24 ) results from shifting X in Eq.( |3.16 ) by an arbitrary (hermitean) matrix while 
Eqs.(CT), ( |3TT91) (or ( gg5[ )) result from the shift X ->• X + e n X n+l . 



3.3 The genus-zero solution 



Brezin and Gross [ |BG80|| with the aid of the Riemann-Hilbert method. The corresponding 
solution had been first found by Kazakov and Kostov [|KK89|| and is discussed in [|GN91 



Substituting this solution into Eq. (|3.15|) which expresses the partition function of the 
Kontsevich model via that for the hermitean matrix in an external field, one gets in genus 
zero 



1 



N 



£ -(A»-2«)*+ii > /A?-2ti + 



u 



A 3 



+ 2N ^ log ^ + A ^' log ( V A ? - 2u + V A ' - 2u 

j 



where u[A] is determined by 



1 \ - 1 

.f(Af-2u) 



(3.25) 



(3.26) 



The solution ( p.25| ) is similar to the strong coupling solution of |BG80|| while Eq. (|3.26j ) 
has emerged to guarantee correct analytic properties. It is important that the r.h.s. of 
Eq. (|3.25|) is stationary w.r.t. u due to Eq. (|3.26|) . For a constant field Aj = (6g)~3, this 
solution recovers the results of Brezin et al. | |BiPZ78 | for the case of a cubic interaction. 

Eq. fl3.26|) can be rewritten in the form of the string equation of a 'general massive 
model' ||BDSS90|1 in genus zero. To this aim, let us expand the r.h.s. of Eq.( |3.26| ) in u 
and substitute 



1 (2n 



N 



nl 







2n-l 



(3.27) 



This equation is nothing but Eq.( [3.9|) with J\f(n) = n\ which fixes the normalization 
|Wit90|| of 2D topological gravity and A = jr as is prescribed by Eq.( p.l4| ). We rewrite 
Eq. (gg|) finally as 



u 



n=0 



(3.28) 



The precise form of the genus-zero string equation can be obtained by the well-known 
shift PW91| : ti->t! + l. 



3.4 Relation to 2D topological and quantum gravities 

It is instructive to compare the solution ( |3.25| ) of the Kontsevich model with known 
results for the partition functions of 2D topological and quantum gravities in genus zero. 
To obtain the perturbative expansion of F [t], one solves Eq. fl3.28p by iterations in u: 

u = t + toh + t tl + t 2 t 2 + ... (3.29) 

with t n ~ tl n+1 , and substitutes the result into the r.h.s. of Eq. ( |3.25| ) which is expanded 



While the complicated structure of the perturbative expansion of Fo[t] represents the 
variety of planar band graphs (taken with appropriate combinatorial coefficients), great 
simplifications occur for derivatives of F [t\. Let us define V by 







(3.30) 



Then, by a direct differentiation of Eq. ([3.25|) , one gets 

VF 



2 ,.2 oo fc+l 

2 _+ 2 h\k + iy 



(3-31) 



k=0 



This expression is again stationary w.r.t. u due to Eq. ( |3.28| ) so that one more application 
of V yields 

V 2 F = u-t (l+t 1 ). (3.32) 

To compare Eq. (|3.25| ) with the known solution of 2D topological gravity in genus zero, 
let us notice that V defined by ( |3.30|) can be rewritten using Eq. ( [3.271 ) as 



V 



71=1 







dt 



(3.33) 



71-1 



This is exactly the operator entering the puncture equation ||DW90|| which reads in genus 
zero: 

dF 



dt 



^ - VF . 
2 



Therefore, one gets from Eq.( |3.31 ) 

dFo 
dt n 



k+l 



u 

y 



(3.34) 



(3.35) 



which is a true formula that gives in particular 

d 2 F Q 



dtl 



u. 



(3.36) 



Since to is the cosmological constant, one sees from this formula u to be the string sus- 
ceptibility. 

Using Eq.( |3.36| ), one can immediately calculate the critical index ^string- For K th 
multi-critical point, when one puts all t n = except for n = and n = K, Eq.( 3.2S ) 
yields u oc xk , ^string = i n full analogy to |[Kaz89|| . Notice, however, that the solution 
( 3.25 ) is associated with the continuum interpolating model while in the standard case 
one 'interpolates' by a matrix model whose couplings should be turned to critical values 
in order to reach the continuum limit. 

Finally, let us mention that the genus-zero solution ( |3.25| ) can be rewritten exactly in 
the form of that for the hermitean one-matrix model in the continuum limit. Let us first 
note that Eq. ( ^-35| ) can be viewed as an integrated version of Eq.( |3.36 ): 



where u(to) is a solution of Eq.( |3.28| ) which is considered as a function of to a t fixed 
values of t n for n > 1. The integration constant is fixed by the fact that u(0) = which is 
nothing but the condition that chooses the perturbative solution of the string equation. 

One more integration of Eq. ( [3.371 ) yields 



Fo 



to 



dx(t Q — x)u(x) + ti 



(3.38) 



which coincides with the representation ( p. 32 ) of the free energy for 2D quantum gravity. 
As has been proven recently |[Wit91| , |MMM91]| , the Kontsevich model obeys the same 
set of Virasoro constraints ( |2.55| ), (|2.56|) as 2D quantum gravity. This demonstrates an 
equivalence of 2D topological and quantum gravities to arbitrary genus. 



4 Concluding remarks 



Loop equations turned out to be a useful tool in studies of matrix models as well as of their 
continuum limit associated with 2D quantum gravity with matter. The point is that loop 
equations are literally the Virasoro constraint imposed on the partition function. In the 
continuum limit, this Virasoro symmetry represents the underlying conformal invariance. 

The appearance of new symmetries of loop equations (as well as the very idea of the 
'double scaling limit' [ BK90 , DS9C , GM90b|| ) is very interesting from the viewpoint of 



multi-dimensional loop equations (see ||Mig83|| ). A step along this line has been done in 



b'KJN 9 1| , |DV V91| , |Goe91| | where the W^-algebras were associated with the continuum limit 
of multi-matrix models. It would be interesting to find an analog of this symmetry for 
multi-matrix models at finite N. 

The fact that the Kontsevich matrix model is a solution of the continuum Virasoro 
constraints (and, therefore, of the continuum loop equations) throws light on the origin 
of Virasoro constraints. It is non-trivial that this matrix model appears as a solution of 



the continuum loop equation. This seems to be analogous to the known property ||Mig83 



of multi-dimensional loop equations which possess solutions differing from the original 
Yang-Mills path integral. 
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